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1 Introduction 



This paper is first from an intended series of papers devoted to problems of photonic band 
gap and the propagation of electromagnetic waves in periodic dielectric structures. The 
question of the existence of a photonic band gap can be rephrased by asking whether 
it is possible to prepare a perfect mirror from a transparent material by drilling (in an 
appropriate way) holes to it. The existence of a photonic band gap has been suggested 
theoretically and recently also been observed experimentally 0. Nevertheless, there 
has been no calculation of band structure by more reliable methods, such as the KKR 
(Kohn-Korringa-Rostocker) or APW (augmented plane wave) methods. Even worse, a 
derivation of the photonic KKR method along lines of has not yet been given . Also, 
both theoretically and numerically, only a limited range of dielectric lattices has so far 
been considered - diamond like lattice, fee and simple cubic lattices In all papers 
referred to, "dielectric atoms" (a compact piece of dielectric in which (and only in which) 
permeabilities are allowed to differ from their host or average values) have been assumed 
to be spherically symmetric. 

In attempting a calculation of the band structure of photons on various dielectric 
lattices with spherical atoms we have tried to generalize the scalar bulk KKR method || 
to the form appropriate for photons |J. In the course of its derivation we have found 
that there persists some confusion in the photonic KKR method. This originates in the 
fact that, unlike the scalar case, in the case of photons fields can change dis continuously 
over discontinuities of electric and magnetic permeabilities. This is the essential difference 
between the Schrodinger and the Maxwell equations. Consequently, derivatives of fields 
may be singular and inward and outward limits of volume integrals over atoms may differ. 
Therefore, one has to carefully distinguish between "inward 1 and "outward 1 formulations, 
i.e., whether, in a given on-shell || surface integral, inward or outward limits of fields 
and their derivatives at the atom boundary are taken. In particular, we shall show that 
the photonic KKR method proposed in does not serve the purpose. In 0, fields in 
integral equations are outward limits of fields with respect to the atom boundary, while 
the scalar product with inward limits of fields is taken. 

These subtleties in derivation of the photonic KKR method have stimulated us to 
write this paper. We hope that this derivation may be both instructive and interesting 
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by itself. A more complete treatment of this subject, together with numerical results on 
photonic band structure, will be given elsewhere |J. 



2 Integral equations for electromagnetic potentials 
2.1 Preliminaries 

We shall consider the stationary macroscopic Maxwell equations in a non-conducting me- 
dium with the current and charge densities equal to zero ||. We shall only consider the 
simplest isotropic case, D(r) = e(r)E(r) and B(r) = /i(r)H(r) where D(r) (B(r)) is the 
electric (magnetic) induction. We shall confine ourselves to monochromatic waves to avoid 
non-local time relation between D(r) and E(r) or B(r) and H(r). We shall also allow for 
complex permeabilities, i.e., for an absorption. In such medium the Maxwell equations 
are symmetric under 

E(r)-H(r), H(r) - -E(r), e(r) - /x(r). (1) 

They can be written as 

V x H(r) +i^e(r)E(r) = 0, (2) 

V x E(r) - % - fi(r) H(r) = 0, (3) 

c 

V ■ D(r) = V ■ B(r) = 0. (4) 

Now, by combining @ and (H) one gets 

V x [V x E(r)] - (c;/c) 2 e(r)/i(r)E(r) 

-^(r)[VKr)]x[VxE(r)] = 0, (5) 

as well as 

V x [V x H(r)] - (cVc) 2 e(r)/i(r)H(r) 

-e _1 (r)[Ve(r)] x [V x H(r)] = 0. (6) 

Due to the symmetry ([I]) of the Maxwell equations, frequently only a single equation 
will be written from a conjugate pair. If /i(r) is set to be a constant, //(r) = fi Q , and 
one parametrizes /x e(r) into its constant uniform value /i e and a spatially varying part 
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v(r) = /i [e(r) — € ], one finds a striking resemblence of (|) and the Schrodinger equation, 
with v(r) playing the role of a potential. Theoretically, this analogy serves as a main 
motivation for the search of a photonic band gap. 

We shall see that to solve the Maxwell equations, even in the case when the potentials 
e(r) and fi(r) can be written as a sum of functions of a single coordinate, is a non-trivial 
task because they do not separate. The reason is that any solution of, say, (^) has to 
satisfy the initial conditions (|j) for the Maxwell equations in a dielectric. Provided that 
e(r) is a non-trivial function of all coordinates then (f|) mixes all components of E(r) 
together, 

V-E(r) = --^[Ve(r)-E(r)]. (7) 
e(r) 

Similarly for H(r) in the case of non-trivial /i(r). 

We shall now turn on to the derivation of integral equations for electromagnetic fields. 
We shall consider either finite or infinite region Q with given boundary conditions imposed 
and with one dielectric atom V s inside it. The generalization to several atoms is trivial. 
Since the group of lattice translations is Abelian, the Bloch theorem extends straight- 
forwardly to vector wave functions. Hence, by setting Q to be a primitive cell and with 
generalized periodic boundary conditions imposed, one can directly turn on to the study 
of electromagnetic wave propagation in dielectric lattices. 

2.2 "Outward" integral equations 

The first attempt to derive integral equations for electromagnetic fields in a dielectric in 
connection with the photonic KKR method is discussed in 0. The derivation is basically 
due to Morse || and starts as follows. Provided r is restricted to the interstitial region, 
r e f2\Vs> the electric (magnetic) field E(r) (B(r)) satisfies the equation 

(V 2 + cx 2 )E(r) = 0, (8) 

where a 2 = u) 2 e fi . Let G CT (r,r') be the Green function of the Helmholtz equation in Q 
with suitable boundary conditions imposed, defined by 

(V 2 + a 2 )G CT (r,r / )=5(r-r'). (9) 
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By using the denning equation (||) for G a (r, r') one is tempting to write the following 
on-shell integral equation in the interstitial region, 

E(r) = / 6(r- r')E(r') dV = 

Jn\v a 




which can be obviously rewritten as the surface integral, 

E(r) = / [E(r') (dS' ■ V')G a (r,r') - G a (r,r') (dS' ■ V')E(r')]. (11) 

Now, assuming Q to be a primitive cell, then due to the periodicity of G a (r, r') E(r'), 
the integral over the cell boundary vanishes and one is left with the surface integral over 
dV s . The same equation also holds for B(r). Thus, formally, both B(r) and E(r) satisfy 
the same integral equation (pTTJ ) outside a dielectric atom V s . What makes the difference 
between them are the matching conditions across a boundary of the atom. 

However, there is one subtle point here. In (|TT|) an outward limit of the electric 
field E(r) and its derivatives is taken. On the other hand, the KKR method requires 
knowledge of either interior or exterior solutions. However, the exterior problem is much 
more involved and cannot be solved usually. Only interior solution can (sometimes) be 
obtained in an explicite form, supposing that there are some symmetries of the problem 
(such as spherical symmetry, for example). Therefore, despite that in the surface integral 
(|ljD the outward limit of the electric field E(r) and its derivatives is taken, usually only 
the inward limit of field and its derivatives is at our disposal. Now, our main task is 
to find integral equations for electromagnetic fields in Q in terms of the inward limit of 
fields and their derivatives at the boundary of V s . Provided that e(r) is discontinuous, the 
discontinuity of the normal component E(r) at the atom boundary is 

n-[E + (r)-E_(r)] = ^E_(r) (12) 

where n is the unit normal vector at the discontinuity. Thus, equations (0-0) inevitably 
contain singular terms. This makes a relation between inward and outward equations 
non-trivial. It is this point that was not taken into account in H. 

It what follows it will be assumed that on discontinuities of e(r) the side limits of of 
e(r) and its derivatives are well defined. The same will be assumed for limiting values of 
fields and their derivatives. 
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2.3 "Inward" integral equations 

In trying to find the "inward' 1 integral equations one can try to establish a direct relation 
between fields and their derivatives on both sides of the atom boundary and rewrite ( |Tlf ) 
in terms of inward limits. It can be shown that like-side limits of e(r) determine relation 
of the side limits of E(r) (0), the side limits of derivatives of e(r) do so for the side limits 
of derivatives of E(r) (see Appendix). Before proceeding this way, the "inward" integral 
equations will be derived by introducing electromagnetic potentials. This way seems to 
be easier and safer in treating singularities that appear. 
If the additive material relations are used 

D(r) = e E(r) +47rP(r), B(r) = H(r) + 47rM(r), (13) 

the Maxwell equations take on the following form 

V x B + iue E = 47rj, 

V x E-iuoB = 0, 
V ■ E = 47rp, 
V-B = 0. 

The current and charge density are defined by the relations 

j = -icuP + cV x M(r), 

p = --v-p, 

where P(r) (M(r)) is the electric (magnetic) polarization, 

P(r) = ^ E(r), M(r) = f^t±H(r). (20) 

The equations (|I4|-|T7|) are formally identical to the equations for the electromagnetic field 
in a vacuum. By introducing a vector potential A(r) and a scalar potential 0(r) such that 

B(r) = V x A(r), (21) 

E(r) = iuA(r) - V0(r), (22) 
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(14) 

(15) 
(16) 
(17) 

(18) 
(19) 



with the Lorentz gauge condition imposed, 

V- A(r) =iue <t>{r), (23) 

equations for A(r) and 0(r) can be decoupled |§. 

In order to simplify our discussion the dielectric medium will be assumed to be mag- 
netically isotropic and from now on we shall set /i(r) = 1 as well as c = 1. For the 
potentials one gets 

V 2 A(r) + cr 2 A(r) = iuo v(r) E(r), (24) 
V 2 0(r) + a 2 0(r) = -V ■ E(r). (25) 
The integral equations for A(r) and <p(r) are 

A(r)=iuf GJr,r')v(r')E(r')d 3 r', (26) 

JVs 

<Kt) = -[ G a (v,v')[V'-E(v')}d 3 v'. (27) 

JVs 

As we have mentioned above, due to the presence of a singular shell at the boundary of 
V s the volume integrals over V s have to be defined carefully as an outward limit through 
measurable sets E n , Q, D S n D V s , lim^oo E n \ V s . The Gauss theorem requires the 
integrand to be at least a continuous function. Provided that a discontinuity appears in 
some region E to which it is going to be applied one proceeds as follows. The region E 
is split (if possible) into non-overlapping subregions Ej, E = UjEj, such that in any Ej 
the hypotheses of the Gauss theorem are satisfied. Then, the Gauss theorem is applied to 
each Ej separately. As mentioned above, inward and outward limits of volume integrals 
with a singular integrand at the boundary will differ in general. First example provides 
Q27|). By using 

V • Kr)E(r)] = -e V • E(r) (28) 

one finds 

0(r):=- / G CT (r,r')[V'-E(r')]dV = 

- / G a (r, r') [V • E(r')] dh' - - I G a (r, r') v(r')(E(r') ■ dS'), (29) 
JVs- ' e Q JdVs- 

where V s+ (K-) means that the outward (inward) limit is taken. Here, the surface integral 

exactly correponds to a delta function singularity which appears due to the discontinuity 
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(|T2|) of the normal component of E(r). Now, it is tranparent that, provided the inward 
limit in fl2~T|) is taken, the resulting integral expressions (|26| - p7D for electromagnetic poten- 
tials do not satisfy the Lorentz condition (p3|). 

The off-shell integral equation for E(r) can be directly written by using the defining 
relations and (|2|0), 

E(r) = -uo 2 I G a (r, r') v(r') E(r') dV- 

-/ [V'a(r,r')][V'-E(r')]dy. (30) 

J v s + 

To find an on-shell integral equation for E(r) we shall make repeated use of @ together 
with the identity (p8j) and with @ satisfied by G CT (r, r'). For a given region E,Ec K, 
one finds in the limit H V s 

u? [ G (r (r,r')«(r / )E(r')d 3 r / -> 

- / [G ff (r, r')(dS' ■ V')E(r') - E(r')(dS' • V^r, r')] 

- / 5(r - r') E(r') dV + / G ff (r,r') V'[V • E(r')l dV. (31) 

Here, the l.h.s. of ( |3"1"D is well defined as a both-side limit. Now, (BUI) can be rewritten 
on-shell as follows, 

E(r) = XVs (r) E(r) + I [G a (r, r') (dS' ■ V')E(r') - E(r')(dS' . V) G a (r, r')] 

[V'G^r, 0^(0 (E^). dS')- 
-l G? <r (^,^ , )[V'•E(0]dS , , (32) 

with Xv s ( r ) the characteristic function of Vg. 

As for B(r), equations ( pT| , p6D together with (fj) imply that 



B(t) 



(33) 



/ ^ [V'G ff (r, r')] x [V x B(r')] dV. 

By repeated use of the Gauss theorem one finds the desired on-shell result, 

B(r) = X v s (r)B(r) + / ^ G a (r, r')dS' x [V x B(r')] 

Jdv 3 - e(r') 

+ / {G CT (r, r') (dS' ■ V')B(r') - B(r') (dS' • V')G CT (r, r')} . 

JdV s - 

We shall now compare outward (JTl|) and inward ([32], integral equations. For 
simplicity, we shall confine ourselves to (52). 



(34) 



2.4 Relation of "out" and "in" formulations 

Here, we would like to derive fl32|) from (|11|) (or vice versa) by direct calculation of relation 
of E(r) and its normal derivative on different sides of E (see Appendix for details). For 
simplicity, we shall confine ourselves to such E for which (d n e) + = (<9 n e)_ = 0. Then the 
last surface integral in (|32|) vanishes. 

As can be found in any textbook, equations (|3|) or (|j) imply the continuity of the nor- 
mal component D n (r) of D(r) (i. e., a discontinuity of E n (r)) or the tangential component 
Et{r) of E(r) on E, respectively, 

K - E- = e -^±E~ = V -^1e~, Et = E~. (35) 

On the other hand, the normal derivative d n E t (r) of the tangential component of E(r) 
changes discontinously across E, 

<9„E+ - d n Et =(—-—) V t D n = -V t E-. (36) 

Rather surprisingly, although the normal component E n (r) is discontinuous, 

d n E+ = d n E~. (37) 

Now, by using (|33| ) one can check directly the equivalence of (O) and (R^) for nor- 
mal components. As for the tangential components, one uses, assuming dV s to have no 
boundary, that 

[V'G CT (r,rOMr')(E(r')-dS') = 

e Q Jdv s - 

-l dS'G a {ry)V'W)E n {r% (38) 
e Jdv s - 

and 

V t [v{r)E n (r)} = v{r)V t E n {r). (39) 



The use of (|36f) and (|37|) then gives the desired equivalence of (0) and (^3), as expected. 



3 Photonic KKR method 



The KKR method is basically the method of rewriting integral equations like p^j3^ ) into 
algebraic ones by using an appropriate basis (a basis of spherical harmonics in our case). 
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In the case of electrons it is known to lead to a very compact scheme if the perturbing 
periodic potential v(r) is spherically symmetric within inscribed spheres, and zero (con- 
stant) elsewhere ||. The band structure of the problem is then determined by purely 
geometrical structure constants, characteristic of the type of lattice under considerations, 
and by phase shifts (logarithmic derivatives, at the surface of the inscribed sphere, of 
the s, p, d . . . radial solutions of the corresponding radial equation), characteristic of the 
scattering properties of a given dielectric sphere. 

In the case of photons the KKR method has been used only within the scalar approx- 
imation to the Maxwell equations |9[]. Recently, a version of the KKR method has been 
given in the framework of multiple- scattering theory [[TO . However, we have obtained a 



different result |pjj| . Nevertheless, a derivation of the photonic KKR method in the spirit 
of is still lacking. 

Let us consider a dielectric lattice and look for Bloch wave-type solutions. In this 
case G a (r, r') possesses the standard expansion in terms of the eigenfunctions of the 
homogeneous boundary value problem, 

C (r r'l - 1 V ex P[*( K n + k )-(r- r/ )] , . 

G ^^)-~-\ (Kn + k)* - <r» ' (40) 

where r is the volume of the primitive cell Q, k is the Bloch momentum, and summation 
runs over all vectors K n of a reciprocal lattice. 

To derive the photonic KKR method we shall turn back to either (|32|) or (0). There is 
one subtlety therein with the characteristic function Xv s ( r )i as wen as i n (FD- Depending 
on the order of the limits taken x(r) may be either one or zero for r e dV s . In order that 
x(r) be zero for r e dV s , one takes as the first limit lim r ^gy s , and then lims^y s _ in (|32j). 
In this case (^) and (^) give an integral equation for E(r) and B(r), respectively. For a 
dielectric lattice of spherical atoms with radius r s it means that the expansion of G a (r, r') 
in spherical harmonics Yi m (6, 0) with r s > r > r' is used, 

Gv(r,r') = H [Aim;iim>ji(o-r')j l ,(crr) + 

lm I'm 1 

+ ffWmm'mH 3i{vr')\ x Y lm (9', <p') Y; m ,(6, 0), (41) 
where ji(x) is the spherical Bessel function defined by 

fa) = (^) Ji+i/2(x), (42) 
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ni(x) := (n/2x) 1 / 2 J-i-i/2(x) , and 9, and 9', (J)' are polar angles of r and r', relative to 
some fixed system of coordinates. This expansion follows by expanding the exponentials 



in ( [TOD for G a {r, r') || |12| . The "structure constants" Aj m; j/ m /, Aj/ TO /.| m = Ai m .^ m *, which 
are functions of cr and k, are characteristic for the lattice under consideration. They are 
exactly the same as for the case of electrons, i.e., as in the case of the Schrodinger equation. 
For r < r s the true solution of our problem, for frequency u, can be expanded in 



vectorial spherical harmonics 12 

oo I 



E(r)=E E -C^Vx [i?f(r)X te (M)] + C,t<(r)X !m ( 



(43) 



i=0 m=-Z 

where the coefficients and specify the amounts of electric (l,m) and magnetic 
(l,m) multipole fields. Rf(r), A = E, M, satisfies the corresponding radial part of 
equation @ or (|B|) 0, with the boundary conditions 

R?(0) = finite, Rf(r s ) = 1. (44) 

Rf 1 and Rf differ in general if e(r) 7^ const. Now one takes scalar products of electric 
and magnetic multipoles with both sides of ( p2|) and obtains a matrix equation. The con- 
dition of solvability, i.e., the vanishing of a corresponding determinant, then determines 
dispersion relation and eventually photonic bands. 

On the other hand, the traditional variational derivation of the KKR method from a 
variational principle uses the order of the limits which makes xv s ( r ) zero f° r r 9V S 
(see (E5t) below). 



3.1 Variational principle for the Maxwell equations 

Let us formulate a general variational principle for the Maxwell equations which holds 
for an arbitrary shape of the basic "atom" of a dielectric lattice. Taking into account the 
off-shell integral equation (|30|) for E(r) the off-shell photonic analogue of the scalar KKR 
functional is defined to be 

A = u 2 I d z vv(v)W(v) ■ (E(r) + u 2 [ GJr, r') v(r') E(r')dV+ 

JVs I JVs 

+ [ V'^^rOtV'-E^dvl. (45) 
Jv a+ J 
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Here, the integral over r is well defined and exists as a bothside limit. Provided e(r) 
is real (no absorbtion), and using well-known hermitian properties of Green functions, 
d T G a (r, r') = — d r iG a (r, r') and G*(r,r') = G^r', r), one can check that variations of A 
with respect to E(r) or E*(r) reproduce correctly the equation for the electric field E(r) 
( P0| ) or its complex conjugate E*(r) within V s , respectively. To deal properly with the 
singularities of G we must use a limiting procedure in evaluating A. In analogy with || 
we set 



A := lim A e , (46) 



where schematically 



A f := u 2 



f d 3 rt>(r)E*(r) • {E(r) + 
+ co 2 f G ff (r,r>(r')E(r') ( lV+ / V'G a {r, r')[V • E(r')] d 3 

+ -<f [V'Gj(r, r')]^(r')(E(r / ) • dS')\ . (47) 

Here, V s - e (dV s - e ) denotes the volume (boundary) of an atom up to a shell of a width of 
e, one side of which is the boundary of the atom and another is formed by the boundary 
shifted by e in the direction of the inward normal. Similarly for V s _2 e and dV s -2 e . 

Supposing that a complete set of interior solutions of (|5j) is known, it is more convenient 



to find the on-shell KKR functional. By using on-shell formulae (|31],[!^) one finds the 
following on-shell form of the variational KKR functional, 

A := lim jjf^ [[(dS ■ V)E*(r)] - E*(r)(dS • V)] - 

/ V[V-E*(r)]d 3 r) ■ (/ UdS' ■ V')E(r') — E(r') (dS' ■ V')] G a (r, r') 

JVs-2, J [JdV s - e 

- - I [V'G^r, r')] v(r') (E(r') • dSf)+ 
-I [G^r.OV'-ECOldS'l. (48) 

JdVs-e J 

Note that unless e(r) is piecewise constant, i.e., the classical muffin-tin potential, it is 
impossible (at least for our variational principle) to write (|45| ) on-shell in terms of surface 
integrals only. Thus, when working with the variational KKR functional it is necessary to 
confine ourselves to the case when e(r) is a real muffin-tin potential. Then the dielectric 
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potential v(r) = v can be moved in front of the integration sign, the terms which contain 
V ■ E(r) vanish, and ([5]) simplifies within V s to 

(V 2 + p 2 )E(r) = 0, (49) 

where p = ui\/e. 

In the case of spherically symmetric dielectric atoms general formula (^) can be 
simplified further, 

A := lim / dS [<9 r E*(r) - E*(r)<9 r ] • 

e-+0 J r =r 3 -2e 




dS' [^E(r') - E(r') d r ] G a (r, r')- 



- - I [V'G a (r, r')] v(r') (E(r') • dS')\ . (50) 

£ Jr'=r s —e J 

The final expression fl50"|) resembles the scalar case |J, the only difference being the last 
term. 

From now on, a further treatment along the lines in is straightforward, albeit 
more involved. The true solution of our problem, for frequency uj and for r < r s , is 
approximated by finite sums of the form ( pffi ) with Rf(r) being now Rf(r) = aiji(pr), 
ai being a normalization constant (see for example |12|]). Then, as in the previous case, 
bands are obtained from the condition of solvability which requires that the determinant 
of A be zero. 

4 Conclusions 

We have presented and discussed the relationship between the outward and inward integral 
equations for electromagnetic waves and we have outlined a derivation of the KKR method. 
For a correct implementation of the photonic KKR method in the spirit of 0, inward 
on-shell integral equations for electromagnetic fields are necessary. The photonic KKR 
method then follows from them either directly or through a variational principle. The 
"direct" photonic KKR method has a wider region of applications : complex and non- 
constant e(r) within atoms. Rather surprisingly, we have found that the variational KKR 
method has a rather limited range of applications. Unless e(r) is piecewise constant 
the variational KKR method cannot be written in terms of surface integrals. Also, for 
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complex e(r), the variational KKR method cannot be used. However, being variational it 
is expected to converge more rapidly within its range of application. Both of the photonic 
KKR methods presented make full use of the same structure constants as the electronic 
KKR method. We hope soon to report on numerical results || . 

In our formulation of the photonic KKR method singularities due to discontinuities of 
permeabilities are safely treated. Our results also imply that the photonic KKR method 
as proposed in can be only used if, at the atom boundary permeabilities are continuous 
and in addition side limits of their derivatives are identically zero. One's first thought 
may be to define the atom boundary to comprise discontinuities of permeabilities such 
that the photonic KKR method of J?| can be used. However, afterwords, singularities of 
delta function are moved into coefficients of corresponding radial equations. 

Multiple-scattering theory for photons, which uses essentially outward integral equa- 
tions, is presented in [ITT] . 
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A Appendix : Derivatives of E(r) on different sides 
of a normal discontinuity 

Let S be some surface discontinuity of e(r) and /i(r) = 1 in O. We shall assume that near 
E only the normal derivative of e(r) may be non-zero. In what follows, such a discontinuity 
will be called normal. In order to determine the discontinuity of <9 n E t (r) we shall look 
carefully at the tangential components of (||), 

V 2 E 4 (r) - ^(r)- 1 V t D n (r) + cA(r)E t (r) = 0, (51) 
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where we have used V • E(r) = D(r) ■ Ve(r). The most singular terms are given by the 
normal derivatives of e(r). Obviously, 



e r 





" i r 


s 





5(r — E) n. 



(52) 



The requirement of cancellation of singular terms gives <9 n Et(r) = e : (r) Vt-D n (r), i.e., 



(0 n E+ - 9 n E t " 



V 4 D n = -v t £- 



(53) 



Since V 4 D n (r) is continuous (see below) one also finds <9 n E^~ = ^ t E~. 

On the other hand, to determine the relation between normal derivatives of the normal 
component of E(r) on different sides of the discontinuity E(r) is replaced by (D(r)/e(r)) 
in (El). After some manipulations one finds (for the normal component) 



V 2 D n + e(r) ( dn^dnj D n + u 2 e{r)D n = 0. 



(54) 



Note that terms proportional V 2 (l/e(r)) have cancelled. Now, the requirement of can- 
cellation of remaining singular terms in (|54]) gives conditions on d n D n (r) and Vt-D ra (r). 
Vt_D n (r) changes continuously across S while d n D n (r) has to be discontinuous. The 
discontinuity is such that d n D n (r)/e(r) change continuously, i.e., 



(55) 



This in turn implies that 

(d n E+ - d n E~ 



-^{d n e) + - -J-(d n e)_ j D n . 



(56) 



Thus, although the normal component E n (r) of E(r) is discontinuous, the side limits of 
d n E n (r) coincide provided the side limits of the derivatives of e(r) at £ are zero. 

Note that the above relations are a general consequence of the (non-stationary) Maxwell 
equations in a dielectric. 
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